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Abstract 



This paper concerns theoretical modeling of degradation of signal with OFDM coding caused by pseudo-random nonlinear distor- 
tions introduced by an analog-to-digital or digital-to-analog converter. A new quantity, effective number of samples, is defined and 
used for derivation of accurate expressions for autocorrelation function and the total power of the distortions. The derivation is 
based on probabilistic model of the signal and its transition probability. It is shown, that for digital (discrete and quantized) signals 
the effective number of samples replaces the total number of samples and is the proper quantity defining their properties. 
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1. Introduction 

+- > ■ 

In the recent years in an enormous number of telecommu- 
nication systems there has been implemented OFDM coding. 
^ This method has been developed already more than 40 years 
. i-H ,ago llrlSD, but for several years has not been used extensively, 
probably partly due to technical limitations. Currently its bene- 
(— i fits are used in systems like wireless computer networks, mobile 
Q_i communication networks, terrestrial digital television or opti- 
L ~ ~~ 1 cal fiber networks. The idea of the method is to encode the 
information on a plurality of subcarriers, which are orthogonal 
J> on the time interval corresponding to a single OFDM symbol. 
This results in immunity to destructive interference in multi- 
path propagation, greatly appreciated in indoor environments 
or cities, but involves specific problems, like e.g. high peak-to- 
. average power ratio and problems with nonlinear distortions. In 
particular, one of nonlinear distortion sources can be digital-to- 
analog (D/A) as well as analog-to-digital (A/D) converters. 
The process of D/A or A/D conversion of OFDM signal ge- 
r* nerally involves some inevitable signal degradation 0-01 • It 
. £h can be analyzed and characterized by means of numerical simu- 
lations, like in J3, IH, but although this way one can obtain 
H precise results for the analyzed cases, obviously this approach 
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does not provide knowledge of the system's behavior in gen- 
eral. For wider understanding of the subject it is much more 
preferred to use an analytic model, showing explicitly relations 
in the system. Some of such models have been already devel- 
oped for signal clipping [9], quantization noise [10] and their 
joint effects 111 U Il2ll . and one of important conclusions from 
them was that degradation due to quantization noise can be re- 
duced by oversampling the signal. However, it seems that little 
is still known about nonlinear distortions in the conversion pro- 
cess, which may be important particularly for systems using 



high order modulations with constellations densely populated 
by symbols. For this reason, the aim of this paper is to provide 
a rigorous, analytic model of pseudo-random nonlinear distor- 
tions based on concepts described in Ill3[ll4ll . These two papers 
have laid foundations for such theoretical model, introducing 
the idea of effective number of samples of a digital (discrete 
and quantized) signal and in this paper this quantity is rigor- 
ously related to degradation of the signal. The derived model 
shows, that like quantization noise, pseudo-random nonlinear 
distortions can be limited by oversampling, but that this limit- 
ing is less efficient. 

The organization of the paper is as follows. The definitions 
and assumptions used for the construction of the model are pre- 
sented in Sec. [2] In Sec.[3]there is briefly discussed the possi- 
bility of decreasing degradation by quantization noise by over- 
sampling of the signal. The assumed decomposition of nonli- 
near distortions into deterministic and pseudo-random parts is 
introduced in Sec. |4] Further, Sec. [5] presents the derivation 
of transition probability for OFDM signal, in Sec. [6] there is 
defined the effective number of samples of the signal and an 
expression for this quantity is derived, and Sec. |7]provides the 
theoretical description of the signal's degradation by pseudo- 
random nonlinear distortions. Finally, Sec. [8] summarizes and 
concludes the paper. 

2. Definitions and assumptions 

The subject of the considerations here is a real, digital signal 
with OFDM coding, converted by a D/A or A/D converter with 
the resolution of n bits. The signal is formed by discretization 
and quantization of an ideal OFDM symbol, given by superpo- 
sition of K modulated subcarriers: 
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x(t) = ^ A* cos(«&f + i 

k=\ 



o. 



(1) 



In general, individual subcarriers represent symbols from arbi- 
trary constellations, defined by amplitudes Ak and phases fa, 
constant during the whole time interval 7s of the OFDM sym- 
bol. For simplicity, it is assumed, that all A* and fa are inde- 
pendent random variables, amplitudes Ak e R+ and have the 
same probability distributions, while each fa satisfies 
/ e 2j0A = o. These assumptions are true for usually used quadra- 
ture modulations, like M-QAM and M-PSK with M > 4 (for 
other constellations, which do not meet the last assumptions, 
like BPSK, just some corrections of numerical factors will be 
needed, still, the presented derivations remain essentially valid). 
The OFDM symbol is assumed to have mean value equal to 
and angular frequencies are particularly chosen as 



a>k — ka>\, with <i>i 



In 
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thus the modulated subcarriers are orthogonal on the time in- 
terval of the OFDM symbol. With these assumptions, the mean 
square {Afj is the same for each k, so that the mean power of 
the OFDM symbol is 



K 



(3) 



The discretized form of the OFDM symbol is the sequence 
of A^s samples 



Xi = x(iT), z = 0,l,...,iVs-l, 



(4) 



with i indexing consecutive samples and T = Ts/Ns being the 
sampling period. The discretization density (oversampling rate) 
is determined by the ratio Ns/K, assumed to satisfy the Nyquist 
criterion: Ns/K > 2. It is known from the central limit the- 
orem, that the signal xt can be well approximated by Gaus- 
sian stochastic process with independent samples, taking value 
x with probability density dependent on the root mean square 
value ©: 

1 (_£_ 



(5) 



Digital devices represent numbers with finite precision. In 
particular, the considered converter with resolution of n bits re- 
quires rounding (quantization) of jc, values to integral numbers 
from the set Z„ = {-2" _1 , . . . , - l}, further called "levels". 
According to the above, the digital representation of the OFDM 
symbol is the sequence of steps 



x q = Xi + A q , with x q e Z„ 



\x] - Xi\ = min |f - x t \ , (6) 



where A q is the quantization error (noise). 

The highest and lowest values of th e sig nal are limited, hence 
clipping of the signal occurs [III [3 • It is convenient to 
relate the power of the signal to the converter's dynamic range 
(clipping level) with the help of the coefficient 



The reasonable practical value is a » 4. Further it is assumed 
here, that clipping is insignificant compared to other degrading 
factors. 

In real converters the levels differ from ideal values, what is 
characterized by so-called differential and integral nonlinearity 
fl. This causes nonlinear distortions of the converted signal. 
If the error for level p is A(p), then the signal with nonlinear 
distortions is 



■A(*?), 



(8) 



where the sequence A^x q ) depends on imperfections of the con- 
verter and also on the converted signal itself. This has consi- 
derable consequences, qualitatively different than quantization. 
This results from the fact, that unlike quantization noise, con- 
secutive samples of the nonlinear distortions A(;t q ) can have the 
same value with nonzero probability. 

3. Quantization noise 

It is known, that (using the introduced notation) if cr » 1, 
then the quantization noise A q is very well approximated by 
white noise with uniform distribution over the interval [— |, ^J, 
i.e. of length the same as level separation [131 . In this case, 
essentially by definition, the whole power of the quantization 
noise 

(( A ?) 2 > = h 

spreads equally over all the A^s samples of spectrum, from 
which only 2K correspond to the signal. Thus, the total power 
of the quantization noise in the signal band is 



cr, 
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Hence, the signal-to-noise ratio for quantization noise: 
cr 2 3 (2"\ 2 N s 



SNRq = — = - 1 — I 



(10) 



(11) 



depends directly on oversampling factor A^s / K, scaling factor a 
and converter's resolution n. In particular, this means, that the 
improvement of quantization noise can be obtained by increa- 
sing the number of signal's samples and that the improvement 
of SNRq is simply proportional to that increase. 

4. Decomposition of nonlinear distortions 

The total error of converter level value can result from: con- 
stant offset and nonunitary gain, which shift and change the 
slope of the converter's transient characteristics, and differen- 
tial and integral nonlinearities 0]. Offset and nonunitary 
gain cause no nonlinear distortions and can be relatively easily 
corrected. These two error sources are further neglected here. 
Then, the remaining error results from nonlinearities and can be 
decomposed into two components: 



(7) 



A(p) = A d (p) + A s (p). 



(12) 



The first component Ad(p) is slowly-varying or deterministic 
part of nonlinear distortions. It can be regarded as systematic 
or regular deviation of converter's transient characteristic, ap- 
proximately constant within adjacent levels, i.e. related mainly 
to integral nonlinearity. The second component A s (p) is the 
quickly varying, pseudo-random or stochastic term, represen- 
ting irregular deviations without any clearly visible pattern. It 
is defined to have vanishing mean value: 



<A,(p)> = 

and no correlation between errors of particular levels: 



(13) 



(A s (p)A s (p')) = A 2 s d pp ,. (14) 

The stochastic part is naturally associated to differential non- 
linearity of A/D or D/A converter. The model derived in this 
paper is taking this source of degradation into account. 

Consecutive samples of pseudo-random nonlinear distortions 
corresponding to digital signal x q form a stochastic process that 
can be written down in the form: 



A.(*?) = J>(p)tf J , 



(15) 
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The autocorrelation function of this stochastic process is not 
trivial, despite of condition (fT4l . In fact, using ([T5T l and then 
(TBI one obtains: 



l? i = (A s (^)A s (^)> = A^(^>. 



(16) 



Obviously for i = the delta is equal to 1 and hence Ro = A 2 . 
For ; + and uncorrected samples with values from a continu- 
ous set, the above delta would determine a zero-measure subset 
of the probability space, hence R, = for i + 0. This is the case 
for quantization noise. However, samples of digital signal be- 
long to discrete set Z„, thus the probability of them being equal 
is greater than zero and autocorrelation function of samples of 
pseudo-random nonlinear distortions is more complicated. Its 
calculation is based on derivation of transition probability for 
the OFDM signal, which is presented in the next section. 

5. Transition probability for digital OFDM signal 

The transition probability of a signal is defined as the proba- 
bility (or probability density for continuous case) of observation 
of given values at two samples with given delay. Exploiting sta- 
tionarity, the transition probability for the digital signal x q can 
be defined as: 



(17) 



For analogue (with continuous values) case, i.e. signal jc, before 
quantization, similarly is defined transition probability density: 

Pi(po,Pi) = P(x =Po A 3d = pi) . (18) 
These two functions are related by the obvious formula: 



Probability of conjunction can be expressed by conditional pro- 
bability, hence: 



r Pi(po,P\) = P(Xi = pi\xo = Po)P(x = p ). 



(20) 



The factor P(xq = po) is given by (0, while the expression for 
P(Xi = P\\xq = po) can be found starting from the definition: 

K 

w(t) = x(t) - xo = J] Re {A*e* (e>*< - l)} . (21) 

k=l 

Then P(x t = pi\ Xo = pa) = P(w(iT) = p x - p \ x = p ). Be- 
cause w(iT) is a linear combination of multiple independent 
random variables, the central limit theorem states, that its pro- 
bability distribution can be well approximated by Gaussian: 



p(yv(iT) = x) = 

with variance 

a 2 wi = (w(jT?) = 2a 2 
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(22) 
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(for i = 1 it will be abbreviated cr w = cr w {). It can be noted, that 
rimj_>o o~wi - 0, just as it should follow from definition of w(t). 
Ignoring any statistical dependence between w(iT) and Xq, the 
sought conditional probability Vi wiiT) = p\ - Pq\xq = po) ~ 
PMjT) = Pl - po), then: 



Pi(P0,Pl) 
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exp 



( Pi ~ P o) 
2cr 2 
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and: 



P?<J>,P')=^— f d A) exp 
27Tcr W iCr J„_i 
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(Pi-Po) 2 ^ 
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This result can be rewritten using the error function 



erf(x) = — — df e '\ 
Jo 



(25) 
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and because cr » 1, it is a good approximation to assume that 
the interval of integration over po is very narrow, the integrand 
is essentially constant within these limits and evaluate it at po = 
p. Then eventually: 



P?(p,p') 
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Thus, the signal value is the same for both samples with proba- 
bility 



P?(p,P r )= [ P+ 'dpo f P+2 dpiP^pi). (19) p* ip>p) *-^ e J-£\ e J—L—\. (28) 



It can be noted, that within the used approximations 



pKp,p) 
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(29) 



and 



-L.^L ? y nxt -_ p) . ,30, 

Assuming, that summation over levels can be replaced by in- 
tegration, the total probability zZp'PoiP' P) ~ 1> meaning, that 
the derived expression is quite accurate. 

Having found the transition probability, the next step is to 
calculate the effective number of samples, defined further in the 
paper. 

6. Effective number of samples of digital OFDM signal 

Because values of xj belong to a set with a finite number 
of elements, the probability that two different samples have the 
same value is greater than zero. This causes, that in the se- 
quence xf appear constant subsequences of two or more sam- 
ples, effectively reducing the number of samples of the signal. 
Let (L) denote the mean length of constant subsequence within 
signal x q . The effective number of samples of the signal x q is 
hereby defined as: 



(31) 



Thus, A^eff denotes the number of changes of values of consecu- 
tive samples in the digital signal. Limiting values are: N e g = A^s 
if no sample is equal to the previous one, and N e s = 1 if the sig- 
nal is constant (all samples are equal). It is convenient to use the 
normalized effective number of samples N e ff/Ns, taking values 
from 1 /Ns up to 1 . 

From the definition it follows, that the effective number of 
samples can be calculated by counting changes of values of 
consecutive samples. Therefore: 



/ JV S -1 , 
JVeff ={l+ 2(1-^)]- 



(32) 



For stationary stochastic process, mean value of the delta in this 
expression is the same for each pair of samples, hence: 

AW 

A^s " V" N S) 



(33) 



The calculated transition probability allows to express the mean 
value of the Kronecker delta: 



!P. q (0,0). (34) 



Then: 
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Again exploiting <x » 1 and approximating the summation over 
levels p by integration, one obtains the closed form expression: 



^,1-erfU- 
Ns \2cr w V2 



\ V2/ V2^2 



(36) 



This formula is compared to results of numerical calculations 
shown in Fig.Q] plotted against resolution n and "reduced reso- 
lution" 

v = n-log 2 — . (37) 

As it can be seen, the derived formula (l36l l reproduces the ef- 
fective number of samples very accurately in most cases. 

The derived formula (l36l can be approximated to show more 
clearly relations between various parameters. For the practi- 
cally significant value of a « 4 one has erf (a/ V5) ~ 1. The 
second term in brackets is inversely proportional to <x, thus 
quickly decreasing as 2~ n . This term is small and can be ig- 
nored. Expanding the trigonometric functions in cr 2 given by 
( T23l into Maclaurin series and then leaving only the leading 
term one obtains 
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For small arguments erf(x) 
simplified formula is: 



2xj y/n, therefore for v > the 
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(39) 



It can be seen that the only dependence on n and Ns/K is 
through the reduced resolution v, what explains the behavior 
of results depicted in Fig.fTbl 

7. Signal degradation by pseudo-random nonlinear distor- 
tions 

The autocorrelation function Rj given by (TTol i contains mean 
value of Kronecker delta, which can be calculated using the 
derived transition probability: 



(40) 



Using (|29l the mean value for samples and i can be expressed 
by mean value for samples and 1 : 



(41) 



Therefore, using (l33l rearranged to express the mean value of 
delta by the effective number of samples, one finds: 



(42) 




Figure 1: Normalized effective number of samples of OFDM signal N e s/Ns vs. fla]) converter's resolution n and (© reduced 
resolution v: comparison of numerical results from NOl and the derived formula (136) . 



Exemplary results of averaged autocorrelation functions from 
numerical simulations and calculated with the help of the above 
formula are presented in Fig. [2] For each case the averaging 
has been conducted over 301 randomly generated signals and 
randomly generated level errors A(p). It can be seen, that the 
derived expression correctly resembles the trends of the consi- 
dered distortions. It is interesting to note, that the autocorrela- 
tion function is generally a peak of certain width, which can be 
characterized by the ratio 



h = i _ 



(43) 



hence, the effective number of samples. It is the expected obser- 
vation, since a lower effective number of samples means more 
repetitions of values in the signal, thus higher correlation bet- 
ween consecutive samples. 

From the Wiener-Khinchin theorem, the power spectrum of 
the pseudo-random distortions can be calculated as Fourier 
transform of Rr. 
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(44) 



The predictions given by (l44l are compared to averaged nume- 
rical results in Fig. [3] A good agreement is observed as the 
theoretical curves resemble the power spectra with high accu- 
racy. The only clearly visible discrepancy is for the frequencies 
with the lowest powers, however it is possible, that the results 
would fit better if the number of results for averaging was in- 
creased. Still, the total contribution of this part of spectrum is 
rather low and even if the derived expression is inaccurate here, 
this will have no significant impact. 
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Figure 2: Normalized autocorrelation function of pseudo- 
random nonlinear distortions: comparison of numerical results 
and the derived formula (l42l for three cases. Values are plotted 
only for a few smallest lags, which are the most interesting part 
of the autocorrelation function. 



The total power of pseudo-random nonlinear distortions in 
the signal's band is 



2nKi 



1 ^ _ Sin A\ 
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(45) 



With R( given by the expression d42t this leads to: 
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(46) 



This expression for cr^. is complicated and it would be bene- 
ficial to derive an approximation that, although less accurate, 
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Figure 3: Normalized power spectrum of pseudo-random non- 
linear distortions: comparison of numerical results and the de- 
rived formula (l44l for three cases. 



would express the relations between quantities more clearly. 
Because Rj is mainly a peak of certain width, it can be approxi- 
mated in (l45l l with a triangular function: 



Rq 
0. 



(R -Ri)i, 



i < 



N s 

otherwise, 



(47) 



where the range of i is shifted by -Ns/2 - this is possible be- 
cause Ri is periodic. Then, nonzero terms of the sum are only 
near i = 0, so one can use sin {2m /TVs) ~ 2ni/Ns and the next 
step is replacement of summation with integration. However, in 
geometrical terms, since in this case summations corresponds 
to adding areas of rectangles, while integration - of triangles, 
and Ro is expected to be the dominant component, to compen- 
sate for ignored area of the central step (peak) a factor of 2 is 
introduced: 2 — * 2 J. Hence, this way one obtains: 
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where Six = x - x i /(3 ■ 3!) + x 5 /(5 ■ 5!) - x 7 /(7 • 7!) + ... is 
the integral sine function. Thus, in the rough approximation 
cr 2 . w /7V e ff and the associated signal-to-noise ratio is 



°~ 1 l 2 "\ N -« 
SNR„ ls = — « — — — . 
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This result has the same form as for quantization noise with the 
effective number of samples /V e ff used in place of total number 
of samples /Vs. 

8. Summary 

The expressions (l42l i and (l44l presented in this paper accu- 
rately describe autocorrelation function and power spectrum of 



distortions caused by pseudo-random, i.e. irregular, nonlinear 
distortions introduced by an A/D or D/A converter into OFDM 
signal and this way extend the theoretical model of the conver- 
sion process present in the literature. Their derivation is based 
on transition probability for OFDM signal and a newly intro- 
duced quantity, the effective number of samples, for which the 
expressions $T7\ and (l36l l have been found. It has been shown, 
that effective number of samples reveals certain information 
about the autocorrelation function and power spectrum of the 
considered distortions. In general, this quantity is important for 
digital signals, because it takes into account important stochas- 
tic property - repetition of values of samples, which has sig- 
nificant consequences. It is expected, that so defined effective 
number of samples should find more applications in description 
of digital signals. 
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